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Abstract. 

The Fluorescence Detector of the Pierre Auger Observatory measures the atmospheric depth, 
O ' Imax, where the longitudinal profile of the high energy air showers reaches its maximum. This 

, is sensitive to the nuclear mass composition of the cosmic rays. Due to its hybrid design, the 

Pierre Auger Observatory also provides independent experimental observables obtained from the 
Surface Detector for the study of the nuclear mass composition. We present X max -distributions 
and an update of the average and RMS values in different energy bins and compare them to 
the predictions for different nuclear masses of the primary particles and hadronic interaction 
^ | models. We also present the results of the composition-sensitive parameters derived from the 

VO . ground level component. 

<N 

1. Introduction 

The determination of the nuclear mass composition of UHECR (Ultra High Energy Cosmic Rays) 
is fundamental to unveil the origin of the most energetic particles known in nature. UHECR 
are detected by means of the extensive air showers created in the Earth's atmosphere, which are 
composed by a cascade of hadrons and electromagnetic (EM) particles. The depth at which the 
EM shower reaches its maximum, X max , strongly correlates with the depth where the primary 
firstly interacted. The X max -distribution carries information about the primary particle and the 
physical processes in the cascade. The hadronic cascade is composed mostly by pions, of which 
the charged pions might decay in flight into a muon and a neutrino. Due to their long lifetime 
and low cross section, muons can leave the core of the hadronic cascade traveling kilometers 
away and be detected. Electrons spread out in time and space because of Coulomb scattering, 
whereas muons practically travel following straight lines. The time structure of the shower disc 
encodes the history of the shower, allowing us to recover information about the longitudinal 
evolution by analyzing the time distribution of the particles arriving at ground level. 

The Pierre Auger Observatory, located on the high plateau of the Pampa Amarilla, is the 
largest cosmic ray observatory ever built. Its hybrid design allows to collect the shower particles 
by a surface detector (SD) and to observe the longitudinal development of the EM profile 
by collecting the UV light with a fluorescence detector (FD). The SD spans 1600 Cherenkov 



detectors in a 1.5 km triangular grid over 3000 km 2 , whereas the FD is composed by 24 telescopes 
distributed over 4 sites overlooking the array. The baseline design is being complemented with 
enhancements both on the SD and FD, and new detection techniques pQ. 

In this paper we present the latest (X max ) and RMS-X max results as a function of energy, and 
the X max -distributions for the highest and lowest energy bins. We also present the evolution with 
energy of SD observables extracted from the time distributions of the signal which are sensitive 
to the longitudinal evolution of the shower. We compare our results with the predictions of 
hadronic interaction models [2]. 

2. Electromagnetic Shower Maximum 

We have used data taken between December 2004 and September 2010, following the analysis 
reported in [3] and updated in [3]. We have considered those showers reconstructed by the FD 
that have at least a signal in one of the SD stations. The longitudinal profile was fitted with 
a Gaisser-Hillas function, and a series of cuts depending on the atmospheric optical conditions 
and on the quality of the reconstructed profile were applied [31 13] • In order to avoid that the 
data sample is biased regarding to the nuclear mass, a number of cuts in the geometry were 
imposed. 

The resolution in X max is around 20 gem -2 . Uncertainties in the atmospheric conditions, 
calibration, event selection and reconstruction result in a systematic uncertainty of < 13 gem -2 . 
The average values as a function of the energy are displayed in the left panel of Fig. [T] (third 
plot from top), alongside with the predictions for different hadronic interaction models. 

The elongation rate Diq = d(X max )/dlog 10 E is better described using two slopes, with 
X 2 /ndf=7A/9, compared to 54/11 if we used a single power law. The change of the elongation 
rate could be interpreted as a change on the nuclear mass composition, provided that there is 
no change on fundamental properties of the hadronic physics. The comparison of (X max ) with 
the different predictions of the hadronic interaction models could be interpreted as a change 
of composition towards heavier elements at the highest energies. The left panel of Fig. [1] also 
displays the RMS-X max , which has been corrected by the detector resolution. It rapidly decreases 
with energy around the same energy of the change on the elongation rate. 

The distributions of X max for all energy bins have been published in [3]. Fig. [2] plots the 
(I max — (X max ))-distributions for the highest and lowest energy bins. The predictions of the 
different hadronic models show a nearly universal shape. At low energies, the shape of the data 
distribution is compatible with a very light or mixed composition, whereas at high energies a 
heavier composition is favored. The high X max tail, characteristic of very light components, is 
suppressed with respect to the lowest energies. 

3. Asymmetry of the signal risetime 

In each SD event, the Cherenkov detectors record the signal as a function of time. The risetime, 
ti/ 2l defined as the time elapsed between the 10% and 50% of the integrated signal, depends 
on the distance to the shower maximum [5], the zenith angle 9 and the distance to the core 
r. In inclined showers, the average risetime depends on an angle defined on the perpendicular 
plane £ as (ti/ 2 /r) = a + 6cos£, due to the different effective distances of the tanks to the 
shower maximum. The evolution of b/a with zenith angle reaches a maximum at O max which is 
different for different primary masses [6J. In the right central panel of Fig. [T]an example of b/a 
as a function of ln(sec#) is shown for the energy bin log 10 (£'/ey) = 18.85 — 19.00. 

Events with energy above 3.16 x 10 18 eV and 9 < 60° were selected jS]. Results of O max as a 
function of energy are displayed in the left panel of Fig. [TJ The systematic uncertainty amounts 
to ~ 10% of the proton-iron separation predicted by models. The number of muons predicted by 
the hadronic interaction models differs from data [9j. A preliminary study including this effect 
indicates a possible change of about < 5% of the proton-iron difference. 
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Figure 1. Left Panel: Results on shower evolution sensitive observables compared with models 
predictions. The error bars correspond to the statistical uncertainty, whereas the systematic 
uncertainty is represented by the shaded bands. Right Panel, top: Typical longitudinal 
development of the muon production. Middle: Average asymmetry in the risetime. Bottom: 
Typical longitudinal development on the energy deposit. 



4. The muon production depth 

An approximated relation between the arrival time delay of a muon with respect to the shower 
front plane and the distance of the muon production point to ground can be used to transform 
the arrival time distributions of muons into production depth distributions [TU]. In [TTj, this 
technique was applied to Auger data in an angular window between 55° and 65° for stations 
with r > 1800 m, as a result of a trade off between EM contamination rejection and the intrinsic 
resolution of the method. In Fig. [TJ right top panel, an example for a real event at E — 94±3 EeV 
is shown. Similarly to the EM component, the depth at which the muon production rate reaches 
a maximum, X^ ax is a good indicator of the primary mass composition. After a series of cuts 
[IT], the measured values of (X^ ax ) are presented in the upper panel of Fig. [TJ The systematic 
uncertainty due to reconstruction bias, core position, rejection of the EM component and quality 
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Figure 2. Centered distributions X max — (X max ), for the lower and highest energy bins. The 
mixed composition corresponds to 50% p and Fe. 



cuts amounts to 11 gem -2 , corresponding to 14% of the proton-iron separation. Note that the 
discrepancy on the total number of muons observed between models and data does not affect 
this result, but only the normalization of the muon production depth distribution. 



5. Conclusions 

Provided that the present models give a fair description of the physical processes and their 
systematics, the average primary composition of UHECR could be inferred from the left panel 
of Fig. [U showing a change towards heavier mass composition at high energies. The evolution 
of (X max ), O max ? and (X^^) is similar in the overlapping energy regions, despite the fact of 
having completely independent systematics. RMS-X max evolution can be accommodated with 
a variety of compositions since it is influenced by the shower-to-shower fluctuations and the 
relative differences on (X max ). The analysis of the shape of the X max -distributions also points 
to the direction of heavier composition as we increase the energy. On the other hand, it is 
not guaranteed that current hadronic interaction models describe well particle physics at these 
energies [9], and a different interpretation would therefore apply. 
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